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Ion Charge States in Halo CMEs: What can we Learn about the 

Explosion? 

Cara E. Rakowski 1,2 , J. Martin Laming 2 & Susan T. Lepri 3 

ABSTRACT 

We describe a new modeling approach to develop a more quantitative under- 
standing of the charge state distributions of the ions of various elements detected 
in situ during halo Coronal Mass Ejection (CME) events by the Advanced Com- 
position Explorer (ACE) satellite. Using a model CME hydrodynamic evolution 
based on observations of CMEs propagating in the plane of the sky and on theo- 
retical models, we integrate time dependent equations for the ionization balance 
of various elements to compare with ACE data. We find that plasma in the 
CME "core" typically requires further heating following filament eruption, with 
thermal energy input similar to the kinetic energy input. This extra heating is 
presumably the result of post eruptive reconnection. Plasma corresponding to 
the CME "cavity" is usually not further ionized, since whether heated or not, 
the low density gives freeze-in close the the Sun. The current analysis is limited 
by ambiguities in the underlying model CME evolution. Such methods are likely 
to reach their full potential when applied to data to be acquired by STEREO 
when at optimum separation. CME evolution observed with one spacecraft may 
be used to interpret CME charge states detected by the other. 

Subject headings: Sun: solar wind — atomic processes — plasmas — waves 
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1. Introduction 

Coronal Mass Ejections (CMEs) represent perhaps the most extreme manifestation of 
solar activity. Of order 10 16 g of plasma is expelled at speeds of several hundred m s^ 1 to 
in excess of one thousand km s -1 . Such speeds are frequently super-Alfvenic with respect to 
the ambient solar wind, and the shocks driven by CME events can be efficient accelerators of 
energetic particles, constituting the main radiation hazard for space-borne instrumentation. 
The spectacular nature of the phenomenon, and its relevance to space based technology, have 
spawned much theoretical and observational work with the ultimate goal of understanding 
CMEs sufficiently deeply to enable the forecasting of such events, in a discipline that has 
become known as "space weather". 

The interplanetary manifestation of CMEs (ICMEs) embedded in the solar wind can 
exhibit a variety of signatures in the magnetic field, solar wind speed and density profiles, 
proton thermal properties, elemental and ionic composition, and the presence of energetic 
particles (see Zurbuchen and Richardson, 2006 for more detail). Of particular interest here 
are the ionic charge states observed in the solar wind during ICMEs. Unlike many plasma 
properties, the charge states are determined within 4 solar radii and remain unchanged as 
they expand further into the heliosphere. This property makes them an important signature 
of the eruption process close to the Sun. A large body of literature exists highlighting 
the rich data sets which detail the unique charge composition observed within ICMEs (e.g. 
Zurbuchen and Richardson 2006, Zurbuchen et al. 2004, Lepri and Zurbuchen 2004, Lepri 
et al. 2001, Henke et al. 2001, Gloeckler et al. 1998, Henke et al. 1998, Galvin & Gloeckler 
1997, Reinard 2005). 

While a number of studies have examined the existing solar wind composition data 
and made inferences on freeze-in temperatures based on computed ionization distributions 
appropriate to coronal equilibrium, in this paper we commence time-dependent modeling 
of the ion charge state distributions of various elements to draw quantitative conclusions 
regarding the thermal energy input and initial conditions in the corona during the CME 
eruption. Using current ionization-recombination calculations and a simple model for the 
spatial and temperature evolution of the CME plasma we attempt to reproduce the charge 
states detected in situ by instruments on the Advanced Composition Explorer (ACE) space- 
craft. In order to provide some context for our work, we first review the status of theory and 
observation for CMEs. Following this, we describe our modeling approach in some detail, 
and apply our methods to a sample of ICMEs detected ACE, which were chosen to provide 
a variety of ICME speeds and charge state distributions. 
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2. CME Theory and Observation: Current Status 



CMEs were first discovered during the Skylab era (jGosling et al.l Il974l ) , and a "halo 
CME " (i.e. an Earth directed event) was first detected by the P78-1 coronagraph ([Howard et al 
19851 ). Since the 1995 launch of SOHO, they are now routinely observed and studied by the 
Extreme Ultraviolet Imagi ng Telescope (EIT ) and Large Angle Spectroscopic Coronagraph 
(LASCO) instruments (e.g. iDere et al.lll997l ). LASCO white light observations are sensitive 
to photospheric radiation Thomson scattered by free electrons in the corona. Three nested 
coronagraphs provided coverage between 1.5 and 30 R@ heliocentric distance, reduced to 2.5 
- 30 R Q following the 1998 demise of CI. EIT records images in narrow EUV wavebands 
emitted closer to the solar surface, and allows investigations of the coronal precursor and 
response to the CME eruption. This instrumentation provides the context for the observa- 
tions of most importance in this work, those of particle charge states and masses detected 
in situ at the LI Lagrange point by mass and charge to mass spectrometers on board the 
Advanced Composition Explorer (ACE). We use Level 2 data supplied by the ACE Science 
Center from the SWICS/SWIMS (Solar Wind Ion Composition Spectrometer/ Solar Wind 
Ions Mass Spectrometer) and SWEPAM (Solar Wind Electron Proton Alpha Monitor) in- 
struments. 

For our purposes, the most important feature of CMEs is their rate of expansion and 
acceleration close to the Sun. For the time being, we are not able to simultaneously observe 
the expansion rate of a CME for which we also detect the particle emission, and must 
estimate the velocity profile of a halo CME from observations of oth er CMEs propagating in 



the plane of the sky. Happily, these frequently follow a similar form (IZhang et al.ll200lL 12004 



Zhang fc Derdl2006l ). Observationally, CME evolution can be divided into three phases. An 
initial phase of expansion at approximately constant velocity in the range 10-100 km s _1 , 
is called the "initiation phase", during which the flux rope rises to a height of about 1.5 
Rq. This is followed by the "acceleration phase". During this period, the CME undergoes 
a roughly constant acceleration up to its final speed, between a few 100 km s _1 to a few 
1000 km s _1 for the fastest CMEs. Typically this final speed is achieved around 5-10 R Q . 
Lastly, the "propagation phase" with essentially constant velocity transports the CME to 1 
AU and beyond. These features are also apparent in theoretical work. In the breakout m odel 
jAntiochoslll998l : kntiochos et alill999l : beVore & Antiochosll200ol ; kulanier et aDl2002f ) , the 
initiation phase corresponds to the initial shearing of the magnetic field lines closest to the 
filament channel in an overlying arcade of loops. As the filament rises and reconnection 
commences under neath it, the acceleration phase begins, typically at 1.5 R & (see Fig. 1 in 
Lynch et al.l l2004f ) . In the breakout model, posteruptive reconnection beneath the erupting 
filament is not strictly necessary for an explosion to occur, though reconnection above the 
filament is required. There appears to be no clear event to signal the end of the acceleration 
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phase, tho ugh clearly a CME c annot continue accelerating forever. Analytic models of CME 
eruptions (jLin fc Forbes! l2000l ) can also give similar height-time and velocity-time profiles. 
Again, the transition from the initiation phase to the acceleration phase corresponds to the 
formation of a current sheet below the flux rope. The resulting analytic solutions for the 
height and velocity of the reconnection-driven flux rope do indeed show an approximately 
constant acceleration phase, followed by a roughly constant velocity phase, at least for models 
with relatively high reconnection rates producing fast CMEs. Some deceleration may also 
occur during the propagation phase. 

The morphology of the erupting CME plasma, while generally quite variable, shows a 
regular pattern of features for flux rope or magnetic cloud CMEs, which are thought to be 
associated with erupti ng prominenc es, and which will be the focus of most of our modeling 
using the event list of iLynchl (120061 ). The CME front (a shock front in CMEs fast enough 
to drive a shock) encloses a "cavity" region , presumably a region of strong magnetic field 
comprising the flux rope (ILynch et al.ll2004l ). Towards the rear of the erupting plasma is 
the CME "core", and region with density perhaps a factor of 10 higher than in the cavity. 
This region possibly has an origin as the cold prominence material. Reconnection above 
the erupting filament in the breakout model would be expected to heat the cavity plasma, if 
anything. Unfortunately the cavity plasma is usually too rarefied for any he ating to be visible 



in th e detected ion charge states. Reconnection behind the eruption (e.g. iRiley et al.l 12002 
20071 ). for example to form postflare loops, will most likely heat the core mat erial, where the 
densit y is high enough that increased ionization can set in before freeze-in. iKumar &: Rust 
f|l996h model a lower velocity flux rope eruption, making the assumption that magnetic 
helicity is conserved during the process. They find that the magnetic energy of the flux 
rope must decrease, and that not all of this energy can be converted into kinetic energy of 
the expanding plasma, some must go to heat, though the exact form and location of the 
magnetic energy dissipation is not specified. 



The flux rope may also be distorted during transit to 1 AU. IRiley &: Crookerl (120041 ) 
model the effects of spherical expansion and the effect of pressure gradients between the CME 
plasma and the ambient solar wind. The flux rope may expand substantially in latitude. 
Schmidt fc Cargilll (120011 ) consider the case of the flux rope propagating in the velocity shear 
layer between the fast and slow solar wind, and find that the magnetic topology may be 
changed by reconnection. In this way, plasma from the CME "core" and "cavity" regions 
may become mixed together before in situ detection at 1 AU by ACE. 
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3. Modeling Approach 

Our simulations follow the evolution of high density core and lower density cavity ma- 
terial for a CME. The ionization balance for the two components are followed separately 
starting from coronal or flare temperatures allowing for heating during the acceleration 
phase of the CME evolution, particularly of the core material. Various parameters listed 
below are adjusted to match the velocity, plasma density, and most importantly the charge 
state distributions of O, Si and Fe observed by ACE for a small selection of CMEs. Given 
the degeneracies in parameter space that will be removed by the simultaneous observation of 
CMEs from face-on and edge-on viewpoints with STEREO, our primary intent at this time 
is to illustrate that the CME charge states can be successfully modeled under reasonable 
assumptions. 



3.1. Charge State Evolution 



Our approach to modeling CME charge states is to follow the behavior of the ionization 
balance of a Lagrangian plasma packet, using an analytic prescription for the hydrody- 
namic or magnetohydrodynamic evolution. We use an adapt ation of the BLASPHEMER 
(BLASt Propagation in Highly EMitt ing EnviRonment)0 code Laming fe GrunlboO^ . 12003 



Laming fe Hwang! 120031 ; lLamingi 12004 ). which follows the time dependent ionization balance 
and temperatures of a Lagrangian plasma parcel as it expands in the solar wind. The ini- 
tial conditions are set by assuming ionization equilibrium at an electron temperature of 
1 — 3 x 10 6 K appropriate for coronal plasma. 



The density rii q of ions of element % with charge q is given by 



dn 



dt 



ion,q^iq 



l+Ti e \C T r,q+\ + Cdr,q+l) n i q+l~ n e {C rTy q + Cdr,q) 



72 



where Cj onig , C rTA , Cd r , q are the rates for electron impact ionization, radiative recombination 
and dielectronic recombination respectively, out of the charge state q. These rates are th e 



same as those used in the recent ionization balance calculations of iMazzotta et al.l (119981 ). 
using subroutines kindly supplied by Dr P. Mazzotta (private communication 2000), with the 
following u pdates. Dielectronic re c ombin ation from H- to He-like and from He- to Li-like are 

taken from iDasgupta fe Whitneyl (12004 ). Dielectronic recombinatio n for the successive iso- 

electronic sequences Li-, Be -, B-, C-, N-, Q-, and F-like are taken fromlColgan. Pindzola. fe Badnell 
(l2004h . lColgan et al.M2003h . lAltun etldl (bo04h . lzatsarinnv et~aH (l2004ah . lMitnik fe Badnell 



lr The name gives away its origin in modeling laboratory and astrophysical shock waves. 
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( 120041 ). IZatsarinny et al.l (120031 ). and iGul (120031 ) respectively. Additionall y, dielectronic re- 



combin ation from Ne- to Na-like and from Na- to Mg-like are taken from IZatsarinny et al. 



from 



Badnell 



(2004b) and IGul (120041 ) . We also take dielectronic recombination rates for Fe 13 
(120061 ). The electron density n e is determined from the condition that the plasma be elec- 
trically neutral. The ion and electron temperatures, T iq and T e are coupled by Coulomb 
collisions by 

(T iq - T e ) q 3 n iq / (q + 1 



■i 'i 



dt 
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dQ 



3n e kB dt 



(3) 



Here M iq is the atomic mass of the ions of element i and charge q in the plasma, and 
In A ~ 28 is the Coulomb logarithm. The term in dQ/dT represents plasma energy losses due 
to ionization and radiation. The term — (T e /n e ) (dn e /dt) ion gives the reduction in electron 
temperature when the electron density increases due to ionization. Recombinations, which 
reduce the electron density, do not result in an increase in the electron temperature in low 
density plasmas, since the energy of the recombined electron is radiated away (in either 
radiative or dielectronic recombination). 



3.2. Hydrodynamic Evolution 



We based our model CME evolution on the phenomenology described in section 2, of 



an initiatio n, acceleration and propagation phase, to use the terminology of (jZhang et al. 



200ll . l2004f ) . Plasma starts at 1.5R & moving at an initial velocity, Vi, between 10 and 100 km 
s" 1 , i.e. as it moves from the "initiation" to the "acceleration" phase. The initial electron 
density, is taken to be either near 10 8 cm -3 or a factor of 10 lower, corresponding to plasma 
in the CME "core" or "cavity" respectively. An acceleration, a, of 0.1 — 0.5 km s -2 is chosen 
so that the CME reaches its final coasting velocity, Vf, at a heliocentric distance of 3 — lOi? 
(i? c ). The density is assumed to fall off as ]_/r( 2+VA '( VA+Vr >> , where va is the coronal Alfven 
speed and v r is the CME expansion speed. This form, with a suitable choice for va (~ 1000 
km s _1 ), gives an initial superradial expansion as the CME expands laterally, going over 
to a 1/r 2 form at large distances from the Sun. The Alfven speed, va, is assumed to vary 
approximately as 1/r 1 / 3 or n l J & , coming from the approximate relations B oc 1/r 4 / 3 (the 
geometric mean of B r oc 1/r 2 and Bq oc B^ oc 1/r) and p oc 1/r 2 . iMann et al.l (120031 ) give a 
more detailed account of the variation of the Alfven speed with radius in the quiescent solar 
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corona above an active region. However our simpler form gives a densit y dependence on 



radius which matches very well with the various plots in lLynch et al.l (120041 ). and is certainly 
adequate for our needs. 

For simplicity we assume the core and cavity differ in density by a factor of 10, but 
contribute roughly equal amounts of material by mass and undergo the same velocity and 
expansion evolution. Their electron temperatures at "initiation" are either assumed to be 
the same or the cavity temperature is held lower, in the range of typical coronal, rather 
than flare, temperatures. Subsequent heating of the higher density core plasma during the 
acceleration phase is explored assuming a heating rate for the CME plasma proportional to 
the rate of kinetic energy increase, i.e. a constant fraction QE/KE during the acceleration 
up to Vf. Here we assume that all particle species are heated equally, but our models are only 
sensitive to the electron temperature. Hence depending on how the magnetic reconnection 
partitions energy between electrons and ions, and whether the heat input is constant during 
the acceleration phase, different CME energy budgets may result. Having the heating be 
proportional to the acceleration during the acceleration phase is one choice which allows for 
particularly easy comparisons of the energetics. However, we cannot distinguish between this 
scenario and impulsive heating to high temperatures (2 — 3xl0 7 K) during the initiation phase 
(with insufficient time to reach ionization equilibrium). The observation of high charge states 
simply requires that high temperatures be reached while the density is still high enough to 
allow for significant ionization of the high Z elements. Such heating was not further explored 
for the cavity material simply because at those densities it would have little or no impact on 
the ionization state which freezes in well below 2R Q . 



4. Modeling a Sample of ACE Coronal Mass Ejections 



A samp le of 6 ICME events was selected from a catalog of magnetic cloud ev e nts b y 



Lynch! (120061 ) , supplemented by two more events from the survey of lUgarte-Urra et al.l (120071 ) , 
all listed in Table 1. Listed here are the average, standard deviation and maxima of the He 
velocity, proton density and the abundance ratios of He/O and Fe/O as measured at ACE 
with the SWICS/SWIMS and SWEPAM instruments (proton density). The abundance 
ratios seen were generally typical for coronal material except for the Halloween 2003 event 
which shows significantly enhanced Fe/O. He/O ratios on the other hand, especially in the 
faster CMEs, are more typ ical of the chromospheric value of ~ 130, or of that found in 
flares (IFeldman et al.ll2005f). and not the lower values foun d elsewhere in the solar wind and 



corona (e.g. 



Laming Sz Feldmanl 120031 ; iKasper et al.l 120071 ). The observed velocity may be 
considered a lower limit, since deceleration is likely to have occurred in transit to 1 AU, 
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however our simulation results were relatively insensitive to the final velocity. The proton 
density at 1 AU can be compared to the electron density at 10 Rq by multiplying by a factor 
of ~ 400 for the 1/r 2 density fall off and another factor of 1.2 — 1.6 for the proton to electron 
density ratio. We expect the measured proton density to reflect some combination of core 
and cavity material. 

Examples of successful CME models for the chosen events are listed in Table 2 in order 
of increasing velocity. The dominant charge states for each model are listed in order of 
abundance composing a total charge fraction of at least 0.8. These are representative of the 
charge states observed at ACE, which, however, did vary over the course of an ICME event. 
Examples of the observed charge states at different slices in time, the final modeled charge 
states of the combined core and cavity plasma, and the initial evolution of ion fractions in 
the core material are shown for the 2001 doy 351 ICME and the Halloween 2003 events in 
figures 1, 2, and 3 and 4, 5, and 6, respectively. The model solution for any given event is not 
unique, and we used this freedom to choose an ensemble that exhibits the range of behavior 
possible and the effect of any given parameter. To begin modeling, the parameters of the 
initial velocity, density, acceleration and lateral expansion (via the Alfven speed) are chosen 
to match the velocity and plasma density of the CME in question as measur ed by ACE. The 



model acce leration and expansion also agree well with observational results in lVourlidas et al. 



( 120031 ) and iThernisien et all (120061 ) for the region behind the forward shock. Both theoretical 
models and observations of CMEs in the plane of the sky show a high density core of plasma 
surrounded by a lower density cavity. The bimodal charge distribution of both Fe and Si in 
most of the CMEs considered here is naturally concordant with a mixing of cavity and core 
material during the CME flux rope passage to 1 AU after freeze-in has occurred. 

In 5 out of 8 sample CMEs the dominant Fe charge states are Neon-like (16+) or higher, 
indicating that high temperatures, comparable to flares (~ 10 7 K), are involved. Starting the 
plasma from this temperature and allowing the ions to recombine as they expand can often 
account for the Fe ionization balance, with peaks around Fe 16+ and Fe 8+ (the Ne-like and Ar- 
like charge states, which have small recombination rates to the next charge states down, hence 
population "bottlenecks" here). However the lower-Z elements place a limit on the maximum 
starting temperature (at least if assuming ionization equilibrium in the seed plasma). Above 
~ 2.5 x 10 6 K O would be mainly 8+ instead of 7+ and 6+ as observed and would not 
recombine significantly during the CME evolution. Evidently plasma must start out much 
cooler, and be further heated as the CME accelerates. While various possibilities exist among 
CME models regarding the role of reconnection in initiating or accelerating the CME, nearly 
all of them requir e reconnectio n to produce the po s t-erup tive arcades of magnetic loops (e.g. 



Riley et al.l 120071 ). We follow iLepri &: Zurbuchenl (120041 ) and argue that this reconnection 



must also heat the CME plasma and that this is the cause for the high charge states. Of 
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the 3 events without a peak at high charge states in Fe, in ICMEs 2003 doy 129 and 2001 
doy 351 the broad distribution of both Si and Fe does require some heating. In the slowest 
event that shows no elevated charge states, 2002 doy 173, no heating is required but neither 
can it be excluded. 

Finally, in Table 2 we also enter for some CMEs the velocity of the CME front as deduced 
from SOHO/LASCO height-time observation^ and the acceleration necessary to produce 
this velocity at the heights observed. The variation in heating rate as a multiple of the CME 
kinetic energy is surprisingly small. Only for the extreme 2003 October 29 (doy 302) event is 
the modeling really sensitive to this behavior. We caution that the height-time observations 
give the velocity of the CME front, not the expansion velocity of the flux rope or magnetic 
cloud containing the ejecta, and that it is not possible in all cases to unambiguously identify 
the disk event that gave rise to the observed ICME. 



Discussion and Conclusions 



Our basic result, that the erupting CME plasma mus t be further heated as it expands, 
is not entirely unexpected from theoretical considerations. Kumar fe Rustl (119961 ) model the 
evolution of an expanding flux rope. For a ratio of initial gravitational energy to magnetic 
energy in the flux rope of 0.22, they find that between 0.58 and 0.78 of the initial magnetic 
energy is converted to forms of energy other the gravitational potential energy and kinetic 
energy of the plasma, i.e., is available for plasma heating or the generation of radiation or 
plasma waves. The maximum amount of magnetic energy that can go to plasma kinetic 
energy is 0.38, so at least 1.5 times the kinetic energy in their model must appear as heat. 



Kumar fe Rustl (119961 ) do not specify the precise mechanism by which the magnetic 



energy is converted into other forms, arguing as they do simply from conservation of total 
energy and magnetic helicity. On the other hand, iLin fe Forbes! (120001 ) considered in a little 
more detail the expulsion of a current carrying flux rope, paying more attention to the initial 
quasi-equilibrium configuration and specifying magnetic reconnection as the mechanism of 
releasing magnetic energy. Their treatment of reconnection assumes, somewhat arbitrarily, 
that all magnetic energy destroyed reappears as plasma kinetic energy, and hence the CME 
speed in their model is t herefore an oy erestimate. However they do comment that previous 
numerical simulations by iForbed (1199 ll ) show that "only about half of the released magnetic 
energy is actually converted into the kinetic energy of the flux rope; the other half goes 
into heating, radiation, and the generation of plasma waves". Therefore these authors also 



2 using the CME catalog at |http: / / cdaw.gsfc.nasa.gov / CME Jist 
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suggest that an amount of thermal energy, similar to the kinetic energy of the CME, should 
be generated, and that this should result from magnetic reconnection. However magnetic 
reconnection does not power the CME. It is accelerated by j x B forces within the flux 
rope, and reconnection eliminates magnetic field that holds the flux rope down in the initial 
quasi-equilibrium configuration. 



Magnetic reconnection plays a similar role in the breakout model of C MEs (lAntiochos 



19981 : lAntiochos et al.lll999l : IPeVore & Antiochosll2000l : lAulanier et al.ll2002h . Again, it plays 
the role of destroying magnetic field that otherwise holds the CME plasma in equilibrium. 
Once this magnetic field is eliminated, the CME is powered by the magnetic field below, 
stressed by shear motions either side of the neutral line. Reconnection at this lower site, 
which forms postflare loops, has nothing to do with allowing th e CME to erupt, but may 
still be a heat source for the CME plasma. Hence in both the iLin fc Forbes! (120001 ) flux 



rope model and in the breakout scenario, no strong prediction exists for how much plasma 
heating should be expected as a function of CME kinetic energy, though simple estimates of 
the magnetic energy destroyed in these cases are indeed of the rig ht order of magnitu de to 
supply the amount of heat we see in the charge state distributions. IShiota et al.l (120051 ) give 
a numerical simulation of a flux rope ejection, and stress the role of the slow mode MHD 
shock in providing heat to the flux rope plasma. 

Our estimates of post-eruptive heating are of course sensitive to the electron heating 
only, although we assume that all plasma parti cles are similar l y hea ted. Observations of post 
eruption electron acceleration are reported by iKlassen et all (120051 ) for the 2003 October 28 
event. Type III, II, and IV radio bursts are detected coincident with the soft X-ray rise 
of the flare. Following that, an impulsive electron event in the energy range 0.027-0.182 
MeV and a gradual electron event with energies 0.31-10.4 MeV were detected by WIND- 
3DP and SOHO/COSTEP, with a total duration of about 30 minutes following CME onset, 
by which time the CME front had expanded to about 6R & . The total energy in these 
ac celerated electrons is significantly less than the heating requirements we give in Table 
2 (IMewaldt et al.l 120051 ) . however the timescale following eruption i s commensura t e wit h 
our requirements. The expansion velocity of the CME front given by IKlassen et al.l (120051 ). 
derived from height-time measurements from SOHO/LASCO observations, is also higher 
than the CME velocity we give in Table 2, and requires a significantly larger acceleration 
to achieve this velocity by about 5.8-R Q where it was first observed by SOHO/LASCO. We 
have experimented with different accelerations, and find that for this particular CME, for 
accelerations above 1-2 km s~ 2 , it is impossible to match the observed charge states for any 
post eruption heating rate. We emphasize that we are principally interested in the plasma in 
the CME "core" , and that the expansion of this material might not be well represented by 
the SOHO/LASCO observations of the CME front. Additionally, this was a halo CME, so 
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SOHO/LASCO observes the CME front off to one side of the main eruption. We anticipate 
that such ambiguities will be substantially resolved with the advent of STEREO data. 

Similar conclusions to ours about the thermal energy inp ut to CMEs have be en reached 
from analysis of ultra-violet spectra taken by SOHO/UVCS. lAkmal et al.l (120011 ) studied a 
480 km s _1 CME observed on 1999 A pril 23, and foun d a th ermal energy comparable to 
the bulk kinetic energy of the plasma. ICiaravella et al.l (200ll ) gave similar results for the 
260 km s" 1 1997 Decemb er 12 CME, while 



event, also examined by lUgarte-Urra et al 



Lee et al.l (120071 ) studied the 2001 December 13 
(l2007h. and included in our Table 2. None of 



these appear in the magnetic cloud event list of iLynchl (120061 ) . 



To summarize, our work on interpreting charge state distributions for the ions of various 
elements support previous ideas in the literature that CME plasma continues to be heated 
as the eruption proceeds. A future quantitative study of this may yield novel insights into 
the mechanism(s) of explosion, especially because the charge states are formed by processes 
close to the Sun, and are then transported unchanged to 1 AU. 

This work has been supported by NASA LWS Grant NNH05AA05I (CER & JML) and 
by an NSF SHINE Postdoctoral Fellowship (STL). It made use of the CME catalog generated 
and maintained at the CDAW Data Center by NASA and The Catholic University of America 
in cooperation with the Naval Research Laboratory. SOHO is a project of international 
cooperation between ESA and NASA. This work has also made use of ACE Level 2 SWICS 
and SWEPAM data provided by the ACE Science Center. 
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Table 1: ICME Observations 



Year 


Start 


v He++ (km 


s- 1 ) a 


Ph+ (cm 


" 3 ) 


He/O b 


Fe/O 






(doy) 


ave. 


max 


ave. 


max. 


ave. 


max. 


ave. 


max. 


2002 


173 


416±19 


464 


4.4±1.1 


6.9 


80±10 


92 


0.15±0.03 


0.19 


2000 


210 


442±34 


474 


15.1±7.5 


35.9 


227±245 


1020 


0.34±0.31 


0.90 


2001 


351 


477±20 


500 


3.6±0.8 


5.5 


123±16 


147 


0.06 


0.07 


2000 


178 


504±42 


569 


5.2±2.3 


13.0 


99±38 


171 


0.21±0.04 


0.29 


1998 


268 


640±77 


793 


3.6±2.2 


11.1 


100±72 


214 


0.28±0.13 


0.56 


2003 


129 


706±78 


855 


3.2±2.2 


10.6 


78±27 


114 


0.18±0.04 


0.25 


2000 


262 


718±47 


804 


4.4±3.1 


13.3 


222±64 


335 


0.28±0.17 


0.67 


2003 


302 


993±305 


1700 


3.1±1.9 


9.2 


168±129 


442 


0.67±0.42 


2.33 



a Ranges given are the standard deviation in the values over the ICME event and do not include uncertainties 
in the measurement. 

b Ratio of the number densities of the given elements 
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Table 2: CME Models 



Year Day Input Parameters Output Parameters 







Vf 

j 


a 


V; 


Dp 
re 


T e 

e 


QE 
KE 


r iurt0 




charge states 






km s _1 


km s -2 


km s _1 


10 7 cm" 3 


10 6 K 




10 3 cm" 3 





Si 


Fe 


2002 


173 


425 


0.1 


10 


10 


1 


0.7 


10 


6 


7,8,6 


8,9 


2002 


173 


425 


0.1 


10 


1 


1 





0.9 


6 


7,8,6 


8,9,10 


2000 


210 


500 


0.1 


20 


10 


2 


5 


30 


7,6,8 


12,11 


16,17,15 


2000 


210 


500 


0.1 


20 


1 


1.3 





2.9 


6 


8,9,7 


11,8,12,10,9 


2001 


351 


500-800 


0.1-0.25 


10 


5 


1.8 


6 


7.0 


6,7 


12,11,10 


14,15,13,16 


2001 


351 


500-800 


0.1-0.25 


10 


0.5 


1.2 





0.8 


6 


8,9,7 


11,10,9,8 


2000 


178 


500-850 


0.1-0.15 


10 


10 


1 


9 


15 


6,7 


10,11,12 


16,15,17 


2000 


178 


500-850 


0.1-0.15 


10 


1 


1 





1.5 


6 


7,8,6 


8,9,10 


1998 


268 


750 


0.1 


15 


10 


2.3 


8 


16 


7,8 


12,11 


16,17 


1998 


268 


750 


0.1 


15 


1 


1 





1.8 


6 


7,8 


8,9,10 


2003 


129 


700 


0.1 


10 


10 


1.2 


2.6 


11 


6 


9,10,8 


13,12,14,15 


2003 


129 


700 


0.1 


10 


1 


1.2 





1.1 


6 


8,9,7 


8,11,10,9 


2000 


262 


700-900 


0.1-0.15 


15 


10 


2.5 


3-2.5 


17 


7,8,6 


12,11 


16,15,14 


2000 


262 


700-900 


0.1-0.15 


15 


1 


1.4 





1.7 


6 


9,8,10 


11,12,10,8 


2003 


302 


1300-2500 


0.2-1 


20 


10 


2.3 


4-8 


15 


7,8,6 


12,11 


16,17,15 


2003 


302 


1300-2500 


0.2-1 


20 


1 


1.2 





1.4 


6 


8,9,7 


8,11,10,9 



"Charge states are listed in order of abundance composing a total charge fraction of at least 0.8. 
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Fig. 1. — Charge state distributions for O, Si, and Fe given in terms of the ion fractions 
for the 2001 day 351 ICME. Plotted here are the ten hour average distributions, offset by 
increments of 0.05. 

Fig. 2. — Model charge state distributions for the 2001 day 351 ICME for parameters given 
in Table 2 with weighting of [0.5,0.5] for the "core" and "cavity" components respectively. 
The bottom panel illustrates the temperature and density evolution close to the solar surface. 



-18- 




1.0000 



0.0001 



Si 9 




Si ,: 



Si 8, 



1.5 2.0 2.5 3.0 3.5 4.0 4.5 
heliocentric distonce (solor rodii) 



1-1 — I — 1 — 1 — 1 — I 1 I 


■ — i — ■ — ■ — ■ — ■ — i ■ ■ — ■ — ■ — i — ■ — ■ ■ ■ — i — ■ — ■ — ■ — ■— | 

Fe' 3 -' 5 * 


P\ z^*- r«'6+ r„'2+ 1 















1.5 2.0 2.5 3.0 3.5 4.0 4.5 
heliocentric distonce (solor rodii) 



Fig. 3. — Evolution of charge state distributions for O, Si, and Fe for the 2001 day 351 
"core" for parameters given in Table 2. The upper right panel illustrates the temperature 
and density evolution close to the solar surface. 
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Fig. 4. — Same as figure 3 for ICME 2003 day 302. Increment between ten hour averages is 
0.1 

Fig. 5.— Same as figure 4 for ICME 2003 day 302, weighting [0.8,0.2] for the "core" and 
"cavity" contributions 
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Fig. 6. — Evolution of charge state distributions for O, Si, and Fe for the 2003 day 302 
"core" for parameters given in Table 2. The upper right panel illustrates the temperature 
and density evolution close to the solar surface. 



